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Abstract 

We investigate the effect of B+L— violating anomalous generation 
of massive right-handed neutrinos on their decouphng, when the right- 
handed neutrino mass is considerably greater than the right-handed gauge 
boson masses. Considering normal annihilation channels, the Lee- Weinberg 
type of calculation, in this case, gives an upper bound of about 700 Gev, 
which casts doubt on the existence of such a right-handed neutrino mass 
greater than right-handed gauge boson masses. We examine the possi- 
bility that a consideration of anomalous effects related to the SU{2)r 
gauge group may turn this into a lower bound ~ 10^ Tev. 

PACS number(s) : 18.80.-k, 14.60. St, ll.lO.Wx. 

I. INTRODUCTION 

Neutrino oscillation interpretation of recent observations of solar and at- 
mospheric neutrino fluxes, although presenting some inconsistencies, may be 
taken to strengthen the idea of non-zero neutrino masses. In this situation, in 
addition to the standard model left-handed neutrinos, the existence and masses 
of right-handed neutrinos assume topical interest. 

The contribution of massive neutrinos to the mass-density of the universe 



allows the setting of a lower bound to such a neutrino mass from the usual cos- 
mological constraints on the age and mass-density of the universe [1,2,3]. The 
standard calculations consider a neutrino mass less than gauge boson masses. 
In the present paper, working in a L-R symmetric extension of the standard 
model [4,5,6], we investigate how the nature of the bound is altered when the 
right-handed neutrino mass is greater than gauge boson masses. 

In these L-R symmetric models, the breaking of SU (2)i{ gauge symmetry 
is associated with a critical temperature. This may, typically, be of the order 
of 1-10 Tcv [7,8,9], and right- handed electron neutrino masses « 10 Tev have 
been considered, yielding a left-handed electron neutrino mass « 10~^° Gev, 
by a see-saw mechanism [9]. Now, B-l-L is not conserved in standard electro- 
weak theory due to an anomaly involving the SU(2) gauge group [10], and, at 
temperatures > 1 Tev, B-l-L— violating transitions occur classically, via thermal 
fluctuations, at rates higher than the expansion rate of the universe [11]. So, 
we may expect that similar anomalous generation of right-handed neutrinos 
(in addition to the left-handed ones), via the L-R symmetric gauge group, may 
become important near the S'?7(2)i{— breaking critical temperature. 

Although there is still a lot of fluidity in the matter, particle physics and 
cosmological bounds usually suggest right-handed Zr and Wr boson masses 
with values > 0.5 TeV and 1.6-3.2 TeV, respectively [12]. If, now, right-handed 
neutrinos of mass > 10 TeV come under consideration in the literature, then it 
becomes necessary to investigate whether anomalous effects can, indeed, modify 
significantly the decoupling of right-handed neutrinos with mass greater than 
right-handed gauge boson masses. 

The plan of the paper is as follows. Section I is the Introduction. In 
Section II, the L-R symmetric model is used to evaluate the reduction rate of 
the right-handed neutrinos, in a standard Lee- Weinberg type of calculation. 
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and to observe how the cosmological bound on their mass becomes an upper 
one, when this mass is greater than the right-handed gauge boson masses. 
In Section III, the anomalous rate of reduction of right-handed neutrinos is 
related to the general anomalous rate of B+L— violating transitions, and the 
qualitative effect of the anomalous rate on the previously obtained mass bound 
is estimated, assuming a generic form for the B+L— violating rate arising from 
the anomaly involving the SU (2)/? gauge group. In Section IV, the influence 
of these anomalous effects on the mass bound is studied numerically, using 
numbers obtained by a simple extrapolation from the SU{2)]^ result. 

II.DECOUPLING WITHOUT ANOMALOUS EFFECTS 
A.BoItzmann equation for processes NN FF. 

We wish to set up a Bolzmann equation for the number density of right-handed 
neutrinos, and, from a calculation of the asymptotic number density, estimate 
the contribution of these neutrinos to the mass-density of the universe, and, 
hence, set bounds to the right-handed neutrino mass [1,3]. 

We will simplify matters by neglecting the decay of right-handed neutrinos. 
Then, the normal electro-weak process chiefly responsible for the reduction of 
right-handed neutrinos may be written as NN FF, where F is a quark or a 
lepton, lighter than N. 

We are interested in investigating the situation when the right- handed neu- 
trino mass is considerably greater than the right-handed gauge boson masses. 

To calculate the rate of reduction of right-handed neutrinos we consider 
the L-R symmetric model [5,6]. This model has pairs of fermion doublets /' 
belonging to different representations of SU(2)iXSU(2)i{XU(l)s_i,, like 
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The numbers refer to the quantum numbers T3l,T3r,B — L, respectively. We 
will also write vr = N,ui, = v. 

The fermion gauge-boson interaction Lagrangian is 

Lint = 9{fl^PLfLr^''L^ + fl^PRTRfW'^ + \9'flA^ - 

where Pl,r = 5(1 ± 75), T is the isospin operator and W^^,B^^ are the gauge 
bosons. The neutral currents are set out in the basis 



A" = sinOiW^'sL + Wsr) +\/cos2eB^', 



Zf" = cose W^3L - sinetanO W^r - tanO/cosW B'', 



Z'^=^l^W^,R-taneB^, 
cost) 

where 6 is the Weinberg angle. 

Neglecting Z-Z' mixing, one gets the Z' neutral current Lagrangian [12] 



^NC — 



9 



e/cosW 



cos 



{sin^9j2f''y''[PLT3L - QsinH]f'+ 



f 



cos^e^fj^^iPRTaR - Qsin^0]f)Z'^. (1) 
/' 

Q is the charge operator. 

The charged current Lagrangian consists of terms of the form 

Lccue = j={i>j^,eLWL+h.c.) + -^{Nj^eRWR+h.c). (2) 

Assuming CP-symmetry, and equilibrium conditions for all relevant particles 
except the N neutrinos, the rate of reduction of N neutrinos per unit volume can 
be obtained from the Boltzmann collision integral for the processes NN — > FF 
[13,14]: 

Ta = ^ / dTTNdTTj^dTrFdTrp{2'K)^S'^{pN +Pn-Pf- Pf) 

F 

l^^r(/jv/jv - fNegfNeg)- (3) 



Here, / is the phase space distribution function and /eg is its equihbrium value. 
|A^:f|^ is the spin averaged matrix element squared, with proper symmetry 
factor, for the process NN — !■ FF, assumed, by CP-symmetry, to be the same 
as that for the process FF NN. The measure diTi = gid^p/{{2-K)^2E), gi 
being the degeneracy number. We assume that there is no significant Fermi 
degeneracy, so that 1 — / « 1. 

Because CP-symmetry has been assumed, we further assume that there is 
no N or TV excess, and we can set n = n, as well as /j.^ = = ^f^. We can, 
then, take 

/iVeg = e-^~/^. 

The summation is over quarks and leptons lighter than N. Let us take i/ 
and N to be electron neutrinos. We assume right-handed neutrinos of the other 
two generations to be much more massive than the N neutrinos, so that they 
are not relevant here. 

It is usual to introduce the thermal average of the annihilation cross-section 
times relative velocity 

< a\v\ >= ^ / dnNdn^diTFdTTp (27r)^5^(pjv +Pn — Pf — Pp) 

"■eg p J 

and write (3) in the form [1] 

r„ =< a\v\ > (n^ - nig), 

where n is the number density of the N neutrinos and rieg is its equilibrium 
value. 

Then the Boltzmann equation for the reduction of N neutrinos by these 
processes, in a universe expanding with R/R = H, becomes [14] 

^+3Hn = - <a\v\>{n^-nlg). (4) 
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B. Calculation of < a\v\ > from L-R symmetric model. 

_ — 

The Feynman diagram for NN — > FF, from (1), is given in Fig. 1. For the 

NN ^ ee amplitude, we get, using (2), an additional contribution from the 




We will work at temperatures T < Tqr, the critical temperature corre- 
sponding to the breaking of SU{2)l x SU{2)r x U{1)b-l to SUi2)L x U{1)y- 
We are going to consider N-type neutrinos with a mass M, which is at least an 
order or two of magnitude larger than Mz' [Mz' > 0.5 Tev [12]). At this energy 
scale, we will approximate all quark and the e,/U, r, v masses to zero (mass of 
topwlTS Gev). 

Next, we assume that u and N have Majorana mass eigenstates [15] 

X = iy + z/'=, u} = N + N'', 

where the superscript "c" refers to the charge conjugate field. 

It is usual to consider a bidoublet and two triplet Higgs particles to generate 
Majorana states [6]. In this paper, however, we do not go into the details of 
any specific model of the Higgs sector. While evaluating the matrix element, 
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we have considered N to be purely Majorana, i.e. we have neglected the contri- 
bution of the vector current and doubled that of the axial current by replacing 
(I + 75) with 275 [12]. 

The spin-averaged matrix element squared, with symmetry factor ^7, for 
the first diagram gives, from (1), 

1 

[{p.k){p.k) + {p.k){p.k) - M'^p.p] 



2 \2 ' 



where, 



Now, 



CvF = T3-2Q sin^e, Caf = 00526*. 



q'' = {k + kf = s = AE^cM 
= 4(M2 + kl^) 
»M|„ 

where {Ecm,^cm) is the 4- vector k in the CM frame. So, we approximate 

l/(g2-M|,)2bylM 

We calculate < a\v\ > in two steps. 
First, we calculate 

If = J diTFdTTp{2n)'^S^{k + k-p-p)\Mj^f 

in the CM frame. The result is 

If = \g'/{64ncos^2e)]{Clp + Clp){2/3)P^ 

where (3 is the relative velocity = 2|kcM|/EcM- This "p-wave" term is a sig- 
nature of Majorana neutrino annihilation. 

In the Lee- Weinberg type of decoupling calculation, the N neutrinos may be 
considered to be non-relativistic, as the relevant temperatures are of the order 
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of M. Then, in the comoving "lab" frame, where k makes an angle a with k, 
Ip = [g^/{MT:cos^29)]{C^p + C\p){2/2,){\l^ + k^ - 2|k||k|cosQ)/M2. (5) 
In the second step we do the thermal averaging. Then, 



< a\v\ >F- 



(27r)^' J (27r)^' 



Calculation, in the non-relativistic approximation for the N neutrinos, gives 



In the case of NN — > ee, the effect of the extra diagram can be taken into 
consideration by the usual Fierz rearrangement, which gives, in this case, 

Cvelcos26 {Cve/cos26) + 1, 

CAelcOs2e {CAelcOs2e) + 1. 

Finally, we get 

<a|^;|>=- ^^^-A^y(C^F + (^lf)- (6) 

So, effectively, < cf\v\ >~ as T ~ M. 

C.Mass Bound for Right-handed Neutrinos. 

Introducing x=M/T and Y=n/s, where s is the entropy density, (4) becomes 
(1.66ff*V^')(MVMp0^5:? = -&lfK < M > {Y' - Yl). (7) 



45 ""^ dx ' 45 



X 



or, 



^ = -0.265*^ < ct\v\ > {MMpi/x^){Y^ - F^). (8) 

We take Mpi = 1.22 x 10^'' Gev, and g* g* 100 just below the critical 
temperature, considering iV, W^,Z'(and N^,Nr) to be massive (we have not 
counted Higgs degrees of freedom). 
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Summing over all quarks and leptons, except the three right- handed neutrinos, 
we get, on calculation, 

^{CIf + C\p) = 3.28 (taking sin^e = 0.23). 

F 

Taking g=0.65, < a\v\ >= 0.01/ (M^a;). 

For massive Majorana neutrinos, we get, in the non-relativistic approximation 
[14], 

Feg = 2.89x IQ-^x^e"''. (9) 



From (8), 



rIY 

^ = -(3.16 X 10i7Ma;3)(y2 _ y^^). (10) 



We write ^ = Y -Y^q. 
Then, before decoupling, Yfa Y^q, and A' ~ 0, giving 

A ^ -Mx^Ylq/[2,.lQ X 10"(2yeg + A)] 

Now, we put A = cYeq at decoupling, where c~ 1. According to the 
general numerical analysis of this type of decoupling [14], c(c + 2) = 2 when 
< (j\v\ >~ T. 

At decoupling, when x = » 1, Y^q —Y^q, and 

A(a;d) ^ cY^M = Ma;|/[3.16 x 10^^(c + 2)]. (11) 

This leads to 

Xd = 35.14 - InM - 1.5 /n(35.14 - InM). (12) 

After decoupling, Y » Y^q and Af^Y. 
From (10), we get 

A' = -3.16 X 10^''AV(Ma;^), 
which gives, on integration, at t — > oo, 

Aoo = ^oo = 2M:cy(3.16 x lO"), 



assuming, Y{xci) >> Foo- 

We will take as our cosmological bound [14] : 

Ojv/i^ < 1, where CIn = Pn/Pc = MsoY^o/pc- 
Here, it is assumed that h> 0.4. 

Taking sq = 2970 cm"^ and po = /i^l.SS x IQ-'^^gcm-^, we get 

= 2.80851 X lO^MFoo (13) 
= 3.62 X lO-^M^a;^, (13a) 

where M is to be taken in Gev. 
At the bound, 

3.62 X lO-^M^a;^ = 1. 

Solving this equation and (12) simultaneously, using simple numerical methods, 

we get 

Xd = 23.55, M = 706 Gev. 

Now, if we omit In M in the third term on the R.H.S. of (12), we get, approxi- 
mately, 

Xd = 29.80 - InM. (14) 

If we make this approximation, the error in Xd is less than 5 percent, even if 
M is as large as 10^ Gev. Using (14), we get 

d{nNh'^)/dM = 3.62 x 10"'' x 2M(29.80 - ZnM)(28.80 - InM), 

which is positive for all practical purposes. 

This means that fl^h"^ < 1 fixes an upper bound for M. 

This can be seen transparently if we work in the approximation 

« Y{xd) « 2Y^g{xd), 
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Table 1 

Mass Bound (no anomalous effects) 

M(Gev) Oiv/i^ 
10,000 160 
5,000 42.6 
1,000 1.95 

750 1.12 

706 1.00 

500 0.516 

250 0.136 



taking c « 1 in (11). 

Then, 

(14) shows that f^w/i^ ~ M2(29.80 - lnM)i, and, so, as M increases, ^l^h 
increases, for all practical values of M. 

This conclusion can be verified, numerically, by giving M different values in 
(12) and substituting the resulting xj, in (13a). The results are shown in Table 
1. 

In the usual Lee- Weinberg case, with M « gauge boson masses, one gets 
a lower bound because < a\v\ >^ M^, which leads to ~ 3 M + constant, 
and O.Nh'^ ~ With M » gauge boson masses, < a\v\ >~ , and this 
makes the difference. 

M < 706 Gev is, in effect, incompatible with our assumption of M > 
right-handed gauge boson masses, because, as we have remarked earlier, parti- 
cle physics and cosmological bounds suggest right-handed gauge boson masses 
~ 0.5 — 1 Tev or more. We have to conclude that the assignment of any real- 
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istic mass, greater than right-handed gauge boson masses, to the right-handed 
neutrinos will violate the cosmological bound flj^h'^ < 1. 

However, we have considered only the normal elecrtroweak process NN — > 
FF. In the next section, we consider, in addition, anomalous processes. 

III. INTRODUCTION OF ANOMALOUS EFFECTS. 
A. Anomalous generation of right-handed neutrinos. 

For the standard model, a classical, unstable, timo-indcpcndcnt solution of 
the equations of motion has been identified [16,17]. This sphaleron solution 
corresponds to the barrier between vacua with different topological numbers. 
A sphaleron-mediated transition over the barrier leads to a fermion-number 
violating transition with \AL\ = 3, \AB\ = 3, of the type 

where a, a' are fermion states, differing by 

I AL| = 3, |AB| = 3, and W'^^ are the initial and final SU(2) gauge boson 
configurations, which are essentially classical. (We are neglecting the small 
effect of the U(l) part [17].) 

All colours and families of quarks and leptons will be generated equally, but, 
in any one transition, only one member per doublet will be found. For the rest 
of this paper, we will neglect family mixing and consider anomalous generation 
for a single (the lightest) family. In this case, |AL| = 1,|AB| = 1. a, a' will 
be restricted by the requirement that the sphaleron must be a colour singlet, 
SU(2) singlet and neutral mediator. There are then two relevant amplitudes, 
which we may write formally as 
< Wf'-uudeWy'- > and < Wfl-uddvW'f >. 

All processes with these amplitudes can occur. For example, a may be the 
vacuum and a' may represent uude or uddu. 
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In the L-R symmetric model, we expect, on general grounds [16,11,9], 
anomalous B+L generation above or just below Tcr, from both SU{2)l and 
SU {2)r gauge boson field configurations with non- vanishing topological charge. 

However, the actual construction of the sphaleron solution depends on the 
details of the Higgs multiplet. The SU{2)l sphaleron [17] was worked out with 
a complex doublet. In the L-R symmetric case, the generation of Majorana 
masses at the higher energy scale results from spontaneous symmetry-breaking 
associated with a SU{2)fi triplet scalar field ( in addition to a SU{2)l triplet 
and a bidoublet which develop v.e.v. at the lower energy scale ) [6]. It has 
been shown [9] that the topological condition necessary for the existence of 
a sphaleron solution is fulfilled for a simplified model of SU{2)fi symmetry- 
breaking at the higher energy scale via a triplet complex scalar field. But, the 
construction of an explicit solution has proved very difficult. 

In this situation, one has to assume [9,18] the occurrence of B-l-L —violation 
via sphalcrons for the SU{2)ii gauge group, in addition to B+L— violation for 
SU (2) L at the higher energy scale. Neglecting mixing parameters between left- 
handed and right-handed gauge bosons, we work with a highly simplified model 
in which the give rise to anomalous generation of leptons and baryons from 
left-handed doublets, and the from right-handed doublets. In particular, 
the will generate, anomalously, right-handed N neutrinos. 

First, we want to relate the rate of production of the right-handed N neu- 
trinos per unit volume to the total rate T of AB = 1,AL = 1 anomalous 
transitions per unit volume. We divide AL = +1 processes into four types 
(assuming distinct flavour eigenstates for N and N): 
I: processes with a N in the final state, e.g. 

\W^'j,vac>-^\w'^tuddN>, 
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I: processes with a -/V in the initial state, e.g. 

\W^'j,Nudd>-^\w'^ivac>, 
m: processes with an e~ in the final state, e.g. 

\W;:^Rvac >^ \Wl^^uude- >, 
ro: processes with an e+ in the initial state, e.g. 

\W^'RUude+ >^ \w'^^vac>. 

Therefore, 

I I m m 

where is a sum over all processes of type i. 

Each process has a rate which is determined in an essentially classical way: 
if the thermal fluctuation has sufficient energy to cross the barrier, the process 
will occur. If iuj~ is the frequency of the unstable (sphaleron) mode, a classical 
statistical mechanics calculation gives [19,20] 

Ti = {u^-/n){ImT/T), (15) 

where T is the free energy. Also, 

(/mJP/T)~e-(^°/^), (16) 

where Vq is the barrier height. 

Because of this essentially classical nature, the barrier- crossing rate, at a given 
temperature, under equilibrium conditions, should be of the same order in 
different channels. In other words, we may expect that the rate of AL = 
1, AS = 1 transitions, featuring one member of a lepton doublet, will be of the 
same order as the rate of such transitions, featuring the other member of the 
doublet. As a flrst approximation, we may take 

^r; + ^r^«^r„ + ^r^«lr. (17) 

I I m m 
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The approximation will be bad when the N neutrinos are way out of equilib- 
rium. In a decoupling study, however, one is interested in finding out when the 
species just falls out of equilibrium. 

Let us now interpret I, I formally as Boltzmann collisional processes 

+ j + — > N + a + b+ 

I: N + a + b + — > t+J + 

CPT ensures that for every process of type I, there is a process of type f with 
the same matrix element Mi. Then, we can write, formally, [13,14] 

Ti — J dTrNdTTadTrb...d7ridTrj ...\A4if 

{2wf6^ipN +Pa +Pb--Pi-Pj-)fj^fah- ■ (18) 

We have again assumed that all relevant species are in equilibrium except the 
right-handed neutrinos, and that there is no significant Fermi degeneracy or 
Bose condensation. Also, 

Ti = j dTrjyd7radTr[,...dT:idTrj...\A4if 

{2n)^S\pjf+Pn+Pi...-p-i-Pj-)fNf-ah- ■ (19) 

In these formal expressions, lAdi]"^ is related to the classical probability and is 
not to be interpreted perturbatively. 

As we are interested here in decoupling and not in baryogcncsis, we will 
neglect small CP-asymmetric effects and assume CP-symmetry. Then, we can 
assume, as in Sec. II. A., 

n = n, Ueq = Ueq, also, /« = 6 t ^ f- = e T etc. (20) 

In this case, we can write, from (18) and (19), 

Ti = Ii neq, and = Ii n, (21) 
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where Ii contains the result of the phase space integrations, apart from Ueq and 
n [14,21]. {Ii = Ij due to (20).) It can be interpreted as a thermally averaged 
width in mode 
From (17) and (21), 

^Ii{n + Ueq) = ^r, and 

I 

yr, = — ^^^r. (22) 

^ 2(n + neq) ^ ' 
For AZ/ = — 1 anomalous transitions, we will get, similarly, a process V with a 

N in the final state, and a process V with a N in the initial state, and a similar 

result 

yTT, = —^ -r'. (23) 

^ ' 2(n + neg) ^ ' 

Now, we are neglecting baryon- and lepton- number excess or deficit. In this 
approximation, we can set hn = 0. We can, then, take [20,22] T = V, i.e. 
the rate of AL = 1 transitions the rate of AL = — 1 transitions. 

Hence, the net rate of reduction of N neutrinos by anomalous processes, per 
unit volume, F^v = rate of such processes, per unit volume, with a N in the 
initial state — rate of such processes, per unit volume, with a N in the final 
state 

I' I 

We finally get 

P _ n — rieq ^ 

2(n + ne,) " ^ ' 
As expected, the anomalous rate vanishes if the N neutrinos are in equilibrium. 

Since, we have Majorana mass states, lepton number violating processes 
with AL = 2 can arise from L-violating terms in the Lagrangian. A variety of 
such processes has been considered in the literature [9,18,23]. These theories 
are, broadly, of two types. 

In one type of theory, a massive right-handed Majorana neutrino is allowed 
to decay. As we are neglecting decay, we have not considered such theories. 
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In the other type of theory, the L- violating term is, effectively, of the type 
{mv/v'^)uuHH, where H is a Higgs scalar, is the v mass, and v is the Higgs 
v.e.v. Usually, in these theories, one takes, jra^ = M. If we consider a similar 
term (1/M)NNHH, we will find a cross-section ~ (1/M^). So, such terms will 
not give processes significantly faster than NN — > FF, with M >>gauge boson 
masses, and may be left out in an exercise where the emphasis is on qualitative 
features of the decoupling. 

Assuming, therefore, CP-symmetry and equilibrium conditions for all rel- 
evant particles except right-handed neutrinos, the rate of reduction, per unit 
volume, of the N neutrinos can be written in the form of a Boltzmann equation 

This is our basic equation. 

We have put F = F/j to indicate that we are considering the anomalous rate 
for right-handed neutrinos. 

B. Anomalous effects in decoupling. 

For T < TcL, the critical temperature for the spontaneous breakdown of 
SU{2)l X U{1)y, F = Fl has been calculated [20,24]. For the right-handed 
case, with T < Tcr, the complication of the Higgs sector has obstructed 
a calculation of Tr. However, the very general considerations mentioned in 
Sec.lll.A. imply that F ~ 6'^°/'^. 

Vo can be estimated heuristically, as follows [25,26] . If we assume a sphaleron 
solution with energy Egp, we can put Vq = Egp. But Egp arises from a classical 
solution, i.e. from a limit where many quanta are involved. We can take the 
energy per quantum Mw, and the average number of quanta --^ {1/aw), 
where aw = Then, Egp is a very general estimate, which should 

hold for the right-handed case also, with Mw = Mwr- So, we can write, on 
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general grounds, 

Tr = R{MwR, T)e-^ M„a/(o^wT) ^ ^26) 

where i?, B depend on the precise form of the symmetry-breaking. However, we 
can say that R will have dimension'~(mass)^, and B will be dimensionless and 
of order 1. Also, the whole idea of separating out the exponential is to isolate 
a prefactor R which can be assumed to vary more slowly (in the left-handed 
case, the prefactor varies as powers of the arguments [20]). 

Introducing x and Y, the Boltzmann equation (25) becomes, from (7), (10), 
and (26), 

dY 1 Y - Y„ 



3.16 X 10^^ 



R'{M, a,x) e-s^/(«»'«) 



writing a = (M/Mwr)- We are considering a>l. 
Compactly, we can write 

— = - R(M a Y-Yeg 3.16 x 10^^ ^ _ yi n /2o^| 

where R! and R are obtainable from R. 

As "a" (the quotes are to avoid confusion with the singular a) increases, the 
first term gains importance because of the exponential. Suppose " a" has a value 
for which the first term is predominant near decoupling. Let us characterise 
decoupling by the simple criterion 

A(a;d) = Y,g{xi). (29) 

This is equivalent to taking c=l in (11). 
Before decoupling, Y w 1^,, A' ~ 0, and 

2(2ye, + A)- ^"^"^ 
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At decoupling, when Y^^ = —Yeg, (30) gives 

[6yeg(ccD)/i?]e^"^/("»"^) = 1, 

and, using the value of Y^q from (9), 

e[^-il^<* G{M,a,Xd) = 1, 

where the prefactor G can be assumed to have a slower variation with and 
"a" than the exponential, because R is a prefactor for which this has been 
assumed. Then, assuming the exponential to dominate, we expect, approxi- 
mately, 

( ^)xci = Bxd ~ constant. (31) 

awd 

Now, the sphaleron decay will produce a N neutrino only if the kinematic 
constraint Egp > M is satisfied. As E^p = BM^rn/aw, so, B/aw > a gives 
an upper limit a' on "a" for anomalous effects to occur. For a < a', B > 0, 
and, if "a" is increased, B decreases, so that, Xd increases. 

We approximate Yoo by Y{xd), so that, from (29), 
Voo ~ ^Yeg{xd)- Then, we get, from (13) and (9), 

Ojv/i^ = 1.62332 X 10*^ x} e"^-* aMwR- (32) 

Since, from (31), a;d ~ ^/(l — ^), the exponential will dominate, and we can 
expect that, as "a" increases, Q,Nh? will decrease. This means that CInH^ < 1 
will give a lower bound on "a", and, hence, on M, for a given Mwr, for a< a'. 

If we can actually find values of the parameter "a"= M/Mwr, within the 
range 1 < a < a', for which the anomalous term in (28) predominates, there 
will not be any hindrance from the Lee- Weinberg type of cosmological bound to 
right-handed neutrinos having masses greater than right-handed gauge boson 
masses. 

So, we find that anomalous reduction of right-handed neutrino number may 
have important effects on the decoupling of such neutrinos. 
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Whether these formal expectations will be borne out depends on the actual 
numbers in F/j. Extrapolating the known result for Fi to the right-handed 

case, keeping wide leeway, we will find that numerical results give cause for 
optimism. 

IV. NUMERICAL RESULTS 

We will take the Fi, given in Reference [20]. 

r.= '";g"^ exp|-l^l. ,33) 

Here, the unstable mode is taken Mw, and Egp = 2{Mw/(xw)E. E 
is a number which depends on (A/5'^):1.56 < E < 2.72 for < A < oo(A 
is the 4-Higgs self-coupling constant). We take E = 2. Mw is temperature 
dependent. 

Mw = Mw{m-{T/Tcf]^, (34) 
and Tc = 3.8Mvy(0) [20]. There is an overall constant k ~ 1 [20,24]. We take 

K = l. 

This expression is valid for 2Mw "C T <C 2Mwlaw- However, the range 
of T may be taken to be Mw « T «; Mw/aw [25] 

We extrapolate this rate to get Tr, in a simple way, using the following 
prescription: 

(i) replace Mw by Mwr, 

(ii) write Tcr = zMwr{0), 

and (iii) include an overall factor b. 

z is not known reliably, because, as yet, there isn't sufficient experimental 
data to evaluate the full L-R Lagrangian, including the Higg sector [27]. For 
large z, (34) shows that Mwr ~ MwrIO)- If z is too small, Mwr will become 
imaginary. We will vary z between 2 and 10. The numerical work will show 
that below z=2, the mass is not real, while there is little change above z=10. 
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Whereas the exponential part in (33) will almost certainly be right for 
Tn (apart from the order one quantity E), the prefactor is bound to require 
considerable modification. Considering the prefactor to be a slowly varying 
quantity, whose main function is to set the numerical scale of the essentially 
exponential variation of Ffl with (1/T), we will allow b to vary from 10~^ — 10^, 
i.e. the decoupling will be investigated with anomalous rates for right-handed 
neutrino reduction varying over 6 orders of magnitude around the rate obtained 
by simple substitution of the right-handed W boson mass in the formula for 
the left-handed case. 

We have, then, 



Introducing x and Y in the above expression, the Boltzmann equation (25) 
becomes 

^ = -f(:c){Y'-Y,l)-gi:c){^^J. (35) 

FVom (10), /(.) = . 

^ . 6 1.53 X , a ,9, ? , 118.98a;,, ,0,9,1, 

Tr gives g{x) = {1 - { — fpexpi ^{1 - {—fP], 

a M\YR[0) zx a zx 

where a = M/Mwr{0). 

For X < Xd, this equation simplifies, as in Section II. C, to 

We choose, again, as an approximate criterion for decoupling : 

A{Xd)^Yeg{Xd) =^ Y{Xd)^2Yeg{Xd). 

At decoupling, Y^^ = — Feg- (36), then, leads to the decoupling condition 
3/(x,)re,(x,) + ^^^ = l. (37) 
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Tablell 
Effect of uncertainty in Tor 



z 


X 


A 


2 


31.5 


41 


3 


31.6 


48 


4 


31.7 


50 


5 


31.7 


52 


6 


31.7 


53 


8 


31.8 


54 


10 


31.8 


54 


100 


31.8 


55 



Again, using the approximation 

Yoo^Y{xd)^2Y,q{xd), 

the cosmological bound becomes, from (32), 

1.62332 X 10'^a;|e-^'*aMvKfl(0) < 1. 

At the bound, 

1.62332 X 10S|e-^<«aMvKfl(0) = 1. (38) 

First, we check the effect of z. Taking My/r{G) = 4000 Gev and b=l, we 
solve (37) and (38), numerically, to obtain values of a;^ = X and "a"=A, for 
which f2jv/i^ is just equal to 1. The results, displayed in Table II, show that, 
as z varies in the range 2 < 2; < 10, X varies from 31.5 to 31.8, and A from 41 
to 54. For z=l, Mwr is no longer real. Also, as expected, z=100 gives for X 
and A practically the same values as given by z=10. 

Having seen that the effect of varying z is small, we set z=4 for subsequent 
numerical work. 
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Table III 

Mass Bound (with anomalous effects) 



a 






100 


234.8 


2.5 X 10"^^ 


75 


72.5 


9.8 X 10-18 


60 


43.4 


1.6 X 10"^ 


50.46 


31.7 


1.00 


40 


22.1 


6.8 X 10^ 


25 


19.6 


4.3 X 10^ 


10 


20.4 


8.2 X 10^ 



We next check that the bound obtained is actually a lower bound. We do 
this by varying "a" around the value A. For each assigned value of "a", we 
solve (37) for Xd, and cvahiato fl^h'^ for this Xd from the LHS of (38). The 
results, displayed in Table III, show that as " a" increases through the value A, 
f2jv/i^ falls through 1, from higher to lower values. 

Finally, we vary b from 10^'' to 10'^ The results arc shown in Tabic IV. 
We find that X changes from 31.8 to 31.6, and A changes from 56 to 46. In 
every case, we have verified the nature of the bound, numerically. The results 
(not exhibited) parallel Table III. The bound remains a lower one. If, of course, 
smaller and smaller values of b are considered, eventually the anomalous effects 
term will be swamped by the second term in (25) or (28), and the bound will 
revert to an upper one. However, the numerical work shows that this does not 
happen even for b=10-^. 

It is necessary to verify that the restriction Mwr < T < Mwr/oiw is 
satisfied. For the lower limit, the worst case occurs when Mwr ~ Mwr{^) = 
4000 Gev. Now, T = 4000 ^/X, and the restriction is satisfied if ^ > X. A 
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Table IV 
Effect of overall uncertainty factor 



b 


X 


A 


0.001 


31.8 


56 


0.01 


31.8 


54 


0.1 


31.7 


52 


1 


31.7 


50 


10 


31.7 


49 


100 


31.6 


47 


1000 


31.6 


46 



perusal of Tables III and IV will show that this is, indeed, so, for the parameter 
ranges considered by us. The stronger restriction, with Mwr replaced by 
2MwR, is, however, not obeyed. 

For the upper limit, the worst case occurs when z, and, hence, M\yr is the 
least. Taking z=2, X=31.5, and A=41, from Table II, we find that Mwr/oiw ~ 
90, 000 Gev, while T 5200 Gev. The restriction is obeyed. 

We check the kinematical constraint Esp > M. As Egp = 2{MwR/aw)E, 
we look only at the case when Mwr is the least, viz. z=2. Egp comes out to be 
> 360,000 Gev, in this case, while, even for A=55, M=220,000 Gev, less than 
Esp, as required. 

V. CONCLUSIONS 

Analysing the decoupling of right-handed neutrinos with mass greater than 
right-handed gauge boson masses, using normal electroweak annihilation chan- 
nels, we find that the cosmological bound f2jv/i^ < 1 leads to an upper bound 
on the right-handed neutrino mass M, of about 700 Gev. What this really 
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means is that a right-handed neutrino mass greater than right-handed gauge 
boson masses is unUkely to be allowed cosmologically. 

If we now assume that anomalous B+L— violating processes work at the 
right-handed symmetry-breaking scale by thermal diffusion over a barrier, sep- 
arating states of different B-l-L, in the same way as this happens at the left- 
handed symmetry-breaking scale, then, we find that it is possible to have a 
lower bound for a right-handed neutrino mass greater than right-handed gauge 
boson masses. 

A numerical extrapolation of the anomalous rate from the lower to the 
higher energy scale, allowing a leeway of six orders of magnitude, confirms 
this possibility. Taking Mwr = 4Tev, a lower bound appears for the right- 
handed neutrino mass at about 50 times the Wr boson mass. However, in the 
absence of an explicit calculation of the anomalous rate for the right-handed 
case, the numbers must only be considered as giving qualitative support to 
the idea that, at Tev energy scales, anomalous generation plays an important 
part in decoupling, and may take away cosmological obstacles to the existence 
of right-handed neutrinos with mass greater than right-handed gauge boson 
masses. To obtain reliable bounds, it is necessary to solve the problem of 
constructing explicitly the sphaleron solution for the right-handed case. 
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